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We study spherically symmetric regular and black hole solutions in the Einstein-Skyrme theory 
with a negative cosmological constant. The Skyrme field configuration depends on the value of 
the cosmological constant in a similar manner to effectively varying the gravitational constant. We 
find the maximum value of the cosmological constant above which there exists no solution. The 
properties of the solutions are discussed in comparison with the asymptotically flat solutions. The 
stability is investigated in detail by solving the linearly perturbed equation numerically. We show 
that there exists a critical value of the cosmological constant above which the solution in the branch 
representing unstable configuration in the asymptotically flat spacetime turns to be linearly stable. 
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1. INTRODUCTION 

The Skyrme model is a unified theory of hadrons proposed by Skyrme P|. The model consists of meson fields 
alone which are represented in terms of angular variables to be multi-valued. Associated with this non-linearity, a 
topological soliton solution called skyrmion arises and its topologically conserved charge is interpreted as the number 
of particle sources, baryon number B. Skyrme discussed the particle nature of a i? = 1 spherically symmetric skyrmion 
by imposing the hedgehog ansatz on the pion fields. Witten showed that QCD in the large-A)- limit reduces to an 
effective theory of mesons, and baryons emerge as solitons in this weakly coupled meson theory The detailed 
analysis for the property of the B = 1 skyrmion as a nucleon was performed in Ref. Q u^n quantization of the 
collective coordinate. Axially symmetric B — 2 skyrmions were found and quantised in Refs. The ground state 

of the solution was shown to have the correct quantum numbers of the deuteron. Remarkably, multi-skyrmions with 
B > 2 possess various discrete symmetries analogously to multi-BPS monopoles |B)I3|' 

It has been known that the Einstein-Skyrme (ESl^swtem pos sesses regular and black hole solutions. The spherically 
symmetric black hole solution with Skyrme hair 1^ llfj lll| and self-gravitating skyrmion [m were investigated. 
It was shown that there exist two fundamental branches of the solutions and interestingly one of the branches 
represents stable configuration under linear perturbations [11113. B = 2 Skyrme black hole and regular solutions 
with axisymmetry were constructed in Ref. [iSj. The review of the black hole solutions with Skyrme hair is given in 
Ref. 13 All of these solutions are, however, constructed in the asymptotically flat spacetime. Recently in Ref. 13 
we considered the Einstein-Skyrme system with a negative cosmological constant and found B = 1 asymptotically 
anti-de Sitter (AdS) black hole solutions. In this paper we develop the previous work of Ref. 13 and study B = 1 
spherically symmetric regular and black hole solutions in the asymptotically AdS spacetime. 

In the context of the Einstein-Yang-Mills (EYM) theory, it was shown that regular and black hole solutions are 
unstable for A > 0 iHlillii. but there exist stable solutions for A < 0 yj, (23 ■ We investigate the linear 
stability of our solutions and discuss in detail to see if the presence of the cosmological constant changes the stability 
properties of the solutions as in the EYM theory. 

There has been an increasing interest in the AdS spacetime. Especially the AdS black hole is an interesting object 
from the holographic point of view in the form of AdS/CFT correspondence (313 Brane world cosmology also 
indicates that there was a period when spacetime was AdS with a negative cosmological constant in the early universe 
(for example, see |l9|, |20l ylj, |23 ). The solutions we obtain in this paper provide a semiclassical framework to study 
the interaction of a baryon and gravity or a primordial black hole with a negative cosmological constant. 


2. THE EINSTEIN-SKYRME MODEL 

The Einstein-Skyrme system with a cosmological constant A is defined by the action 

s = I 2A)-h :^g'^'^tr(L^L^)-k y^5^V‘"tr([L^,Lp][L,„L<,])| (1) 

where = U^d^U and U is an SU{2) chiral field. EV is the pion decay constant and e is a dimensionless parameter. 
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We require that the spacetime recovers the AdS solution at infinity and thus parameterize the metric as 

ds^ = + Cir)-^dr^ + r^(de^ + sin^ 6»V) 


( 2 ) 


where 


C{r) = 1 - ^ . 


(3) 


The topology of AdS spacetime is x and hence the timelike curves are closed. This can be, however, unwinded 
if we consider the covering spacetime with topology 

The B — 1 skyrmion can be obtained by imposing the hedgehog ansatz on the chiral field 


U (r) = cos /(r) + in ■ t sin /(r). 


(4) 


Introducing the dimensionless variables 


X = eF^r , fj.{x) = eFTrGm{r) , A = A/e^F^ 


with 


one obtains the skyrmion energy as 

Es = dTT— J 

' ‘‘*T / + 3^”) 




sinV 

2x2 


2Cf 


sin2/ 


e^x^dx 


(5) 

( 6 ) 

(7) 

( 8 ) 


where we have defined m = x2 + 8 sin^ / and v = sin^ f{x^ + 2 sin^ /). The prime denotes the derivative with respect 
to X. The covariant topological current is defined by 


1 

= -^^^=%^{U-^d,UU-^dpUU-^d,U) 


(9) 


whose zeroth component corresponds to the baryon number density 

go- 1 

27r2 r2 

We impose the boundary condition on the profile function as 

/(x) ^0 as r ^ oo , 


( 10 ) 


( 11 ) 


which ensures the total energy ® to be finite. Then the baryon number becomes 

B= [ d^x = -- [ sin^/d/= T(2/o-sin2/o). (12) 

J tt Jf^ 27r 

For regular solutions, /o = /(O) should take tt in order for the baryon number to be one. For black hole solutions, 
fo = f{xh) with event horizon Xh takes the value less than tt, which means that the solution possesses fractional 
baryonic charge. In this case f{xh) is a shooting parameter determined numerically so as to satisfy the desired 
asymptotical behavior in Eq. CU- 

The Einstein equations with a cosmological constant A takes the form 


( 13 ) 
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which reads 


Goo = S^GToo - Agoo ^ l-C-C'x=^ (cuf^ + J 


[e^^cy 


Gii — SttGTii — Agii 


-1 + G + 


x=y-[ Cuf^ - ^ ) - Ax- 


(14) 

(15) 


where we have defined the coupling constant a = AttGF^ . Consequently following two equations are obtained for the 
gravitational fields 


S' = , -(Gx)' + 1 = ^ C Gm/'^ + ^ ) + ■ 


4x 4 \ - X 

Taking variation of the static energy with respect to the profile /(x), one can get the field equation for matter 


(16) 


f" = 


1 


e^Cu 


-{e^Cuyf'+ 4G/'2 + 1 + 


4sinV\ 5 


e° sin 2/ 


Thus the coupled field equations to be solved are given by 


a 




/" = 


e^Cu 


-{e^CuYf + 4G/'2 + 1 + 


4sinV\ 5 


e° sin 2/ 


(17) 

(18) 

(19) 

( 20 ) 


3. BOUNDARY CONDITIONS 


The boundary conditions for regular solutions are determined by expanding the functions around the origin x = 0 
and comparing the coefficients in each order of x in the field equations flB-CT. As a result, we find 


TT + flX + O(X^) 

(21) 

^0 + ^/?(l + 8AV + O(x3) 

(22) 

|/?(1 + 4/2)x3 + 0(x4). 

(23) 


where /i and So are shooting parameters. /i is chosen so as to satisfy Eq. EJ and So is chosen so as to recover the 
AdS spacetime asymptotically, that is, S(x) —> 0 as x —> oo. 

Similarly, in order to determine the boundary conditions on the regular event horizon, let us expand the fields 
around the horizon x = Xh 


f = fh + fi{x - Xh) + 0{{x - Xh)"^) 

S = Sh + Si(x - Xh) + 0((x - Xh)^) 

M ^ - Xh) + 0((x - Xh)^) . 


(24) 

(25) 

(26) 


Inserting them into the field equations iHHji-® and comparing the coefficients in each order of (x — Xh), one obtains 


fi = 


X? + 4 sin^ fh 


Xh{xl + 8sin^ fh){l - 2ni - Ax^) 


sin 2fh 


^1 = + Ssin^ fh)fi 

4x?i 


1^1 = 


1 + 


2 sin^ fh 


sin^ fh : 


(27) 

(28) 
(29) 


where fh and Sh are shooting parameters with the desired asymptotic behavior /(x), S(x) ^ 0 as x ^ oo. 
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4. NUMERICAL RESULTS 

The profile functions of regular solutions are shown in Fig. Q] for several values of |A| with a = 0.1 fixed. There are 
two branches of the solutions for each value of the cosmological constant as well as the coupling constant. We define 
the solution with larger value of /i as an upper branch and smaller value of fi as a lower branch. The skyrmion shrinks 
as |A| becomes larger in the upper branch and expands slightly in the lower branch. Similar behavior is observed for 
the solutions in the asymptotically flat spacetime when a is increased Thus the variation of cosmological constant 
gives a similar effect on the skyrmion as the variation of the coupling constant does. We found the maximum value 
of the cosmological constant with a > 0 above which there exists no solution. When a = 0, regular solutions exist 
for all values of A < 0. The maximum value of |A| is plotted as a function of a in Fig. |21 |A| decreases monotonically 
as a increases and at a ~ 0.162 only asymptotically flat spacetime solution exists. For a > 0.162, we found no 
solution with A < 0. Fig. 0 shows the dependence of the ADM mass M = fJ,{oo) on a and |A|. The ADM mass 
gives the total energy available in the spacetime and therefore, for regular solutions, it is equivalent to the skyrmion 
energy Es = m(oo) where Es = eEslE-j^ in Eq. Q. As to be expected, the mass increases as a and/or |A| increase. 
It is observed that the presence of the cosmological constant affects the upper-branch significantly more than the 
lower-branch. 

For black hole solutions, the proHle function numerically computed for several values of |A| are shown in Fig. 
The horizon radius and the coupling constant are hxed with Xh = 0.1 and a = 0.02. There are two branches of 
solutions for each value of the cosmological constant as was seen in the regular case. The dependence of profiles on 
the cosmological constant is also similar with the regular case. In the lower branch, however, the change in size is 
much smaller and is almost unrecognizable. 

The black hole skyrmion mass-horizon radius relation is shown in Fig. 0 The skyrmion energy Eg is related to the 
black hole skyrmion mass Mhh = /r(oo) by 

Eg = Mbh - ^i{xh) = ^ . (30) 


Since the entropy of the black hole is written by 


^ ^ (xl\ 
4tiG He^ \ a / 


(31) 


one can see that the upper and lower branch correspond to the high- and low-entropy branch respectively. The 
cosmological constant reduces the entropy of the black hole. The reduction of the entropy is also seen when the 
coupling constant increases as is inferred from Eq. m- 

Eig.inishows the parameter fh as a function of |A| for a = 0.0, 0.02, 0.04 with Xh = 0.1 fixed. The value of fh is 
directly related to the baryon number as can be seen from Eq. m- Thus, in the upper branch, the baryon number 
becomes smaller as |A| becomes larger, which represents the baryon more absorbed by the black hole. On the other 
hand, in the lower branch, the baryon number slightly increases as |A| becomes larger. This result also shows that 
the cosmological constant gives a similar effect on the skyrmion as the coupling constant. 

We found the maximum value of |A| above which there exists no black hole solution for each value of the coupling 
constant. In Fig.jTJ the maximum value of | A| is shown as a function of a. The maximum value decreases monotonically 
as a increases, and at a = 0.126 it becomes zero. Thus at a = 0.126, the asymptotically AdS solution does not exist 
and only the asymptotically flat solution exists. For a > 0.126, we found no solution with A < 0. 

Let us denote that for both the regular and black hole cases, the lower and upper branch solution coalesces at the 
maximum value of |A|. 


5. LINEAR STABILITY ANALYSIS 

In this section we shall examine the linear stability of the solutions described in the previous section. Let us consider 
the time-dependent small fluctuation around the static classical solutions /o, (5o and /xq 



= fo{r) + fi{r,t) 

(32) 

6{r,t) 

= (5o(r) -h(5i(r,<) 

(33) 



(34) 
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From the time-dependent Einstein-Skyrme action 


e^C 


f^ + Cnu + v 


e^dx , 


one obtains the time-dependent field equation as 

s 




e^Vf 


e^C 


uf 


where we have defined Uf = Su/Sf and Vf = Sv/6f. The dot denotes the time derivative. 
The time-dependent Einstein equations are then given by 


Goo = SttGToo ^ 
Gn = SttCTii ^ 
which reads the following two equations 


1 - G - G'cc = 


a 
4 

2Sr^\i 


1 


e2«G 


f + Cr]u + 


2v 


g25 4 


1 


X 

f/2 


e2-5G 


f- + Cr\u-^ 


S' = ^ 


1 


f+r^]u 


4:x 

— (CxY -I-1 = — 7 tv + CS'x. 

2x^ 

Substituting Eqs. into Eqs. (ESl and iTO) gives the linearized equations 

= ^(2uo/o/{+M/ 0 / 0 V 1 ) 

-(e‘^°Gia;)' = -^e^°Vfgfi + e^°CQ5[x. 
Eq. 14111 and the classical field equation Eq. Idtill which can be rewritten as 

^ = (e'«Go«o/^)' - , 

are inserted into Eq. (P|) and resultantly one gets 


-(e^«Gix)' = -(e^«Gouo/^/iy. 

This equation can be integrated immediately to obtain 

Gi = -^CoUofUi • 

Similarly let us linearize the field equation Using Eqs. (gu, gni and gHl , one arrives at 

(e'^°GoUo/i)^ ~ Uofi = g g-, uofi 
e^oCo 


where 


Uo = -(e^“Gow/„/'y+(^e^'>Gow2/'2)'-^e^«GoWo/, 


2 ----- 

Setting /i = ^{x)e^‘^^/g/uo , we derive from Eq. (BtHl as 

1 


-(e^«Go?')' + 


2\/'Uo 


e^"G, “0 


y/uE 


+ —Uo 
Uo 


e = a;^ 


e^oGf 


■e- 


(35) 


(36) 


(37) 

(38) 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 


(47) 


(48) 
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Let us introduce the tortoise coordinate x* such that 

dx* _ 1 

dx e^°Co 

with —oo < x* < + 00 . Eq. (P|) is then reduced to the Strum-Liouville equation 

where 


(49) 


(50) 


Uo = e^°Co 


1 

2,/uE 



Up \ 


+ —Uo 

Up 


(51) 


The classical solution is linearly stable if there exists no negative eigenvalue since imaginary uj represent exponentially 
growing modes. Unfortunately the potential Up has a complicated form and we are unable to discuss the stability 
analytically. Thus we solve the wave equation TO numerically under the boundary conditions that ^ vanishes at 
the boundaries, which ensure the norm of the wave function to be finite. The ground state corresponds to the wave 
function with no node. The nth excited state corresponds to the wave function with n nodes. 

We show the eigenvalue of the ground state as a function of |A| in Fig. IHlfor regular solutions and in Fig. Elfor 
black hole solutions. For all values of |A|, upper branch solutions have positive eigenvalues, and hence they are stable. 
Remarkably the eigenvalues of lower branch solutions increase as |A| increases, and at some value, they cross zero 
to become positive. These figures show clearly how the cosmological constant tilts the eigenvalues to the positive 
direction. Thus lower branch solutions change their stability at the critical value of the cosmological constant. The 
stability of black hole solutions more strongly depends on the cosmological constant than that of regular solutions. 

The results we have obtained indicate that the presence of a negative cosmological constant stabilizes both of regular 
and black hole solutions. This is consistent with the study of the EYM system where it was shown that the EYM 
solution is unstable for A > 0 muilii , but the stable EYM solution exists for A < 0 mm. 

Another important feature for the ES solutions with A < 0 is that only discrete modes exist in both of the branches. 
It is because the potential behaves asymptotically as x^, meaning all the eigenvalues are discretized analogous to the 
harmonic oscillator eigenvalues. 


6. CONCLUSIONS 

We have studied regular and black hole solutions in the Einstein-Skyrme system with a negative cosmological 
constant. There exist two fundamental branches of the solutions. For black holes, these corresponds to the high- 
and low-entropy branch respectively. The increase in the absolute value of the cosmological constant |A| gives similar 
effects on the skyrmion as increasing effectively the value of the gravitational constant. The skyrmion shrinks in the 
upper-branch and expands in the lower-branch as |A| increases. Particularly, in the black hole case the baryon number 
is more absorbed by the black hole corresponding to the increase in |A|. There is the maximum value of |A| above 
which no solution exists for each value of the coupling constant. We have observed that the critical value of a for 
black hole solutions to exist is 0.126 and 0.162 for the regular case. 

The linear stability was examined in detail by solving the linear perturbed wave equation numerically. In the 
asymptotically flat case, it was shown that the upper branch is stable and the lower branch is unstable [ulll^- In the 
AdS case, however, there exist stable solutions even in the lower branch depending on the value of |A|. We have shown 
by numerically solving the linearly perturbed equation that the presence of a negative cosmological constant lifts the 
eigenvalues from negative to positive values. The observation that the negative cosmological constant stabilizes the 
ES solutions is consistent with the results in the EYM theory where stable regular and black hole solutions exist only 
when A < 0 mmiii Let us g ive a brief comment on the catastrophe theory applied to the stability analysis of 
non-abelian black holes in Ref. |29|. It seems that for asymptotically non-flat spacetimes, the catastrophe theory is not 
applicable to the stability analysis of solutions. Although concrete analysis should be performed for any statement 
on this matter to be confirmed, we suspect that the cosmological constant needs to be taken into account as an 
additional control parameter in the parameter space, which makes the dimension of the Whitney surface three with 
S = S{M,Bh,A) in the notation of Ref. |^. This extended catastrophe theory may be applicable for spacetimes 
with A. 
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The solutions exhibited in this paper provide a semiclassical framework to study the interaction of a baryon and 
gravity or a primordial black hole in the presence of a negative cosmological constant. If the universe had gone through 
the AdS phase in the early epoch as indicated in Refs. those solutions may have been produced after 

the hadronization. Since our model predicts the baryon decay by primordial black holes, the Einstein-Skyrme theory 
should be useful as a simple framework to study such decay process. 

The Skyrme model corresponds to QCD in the large N limit and therefore it may be worth understanding the 
Skyrme model in the context of string and brane theories. Also interesting applications are to find solutions with 
nonspherical event horizon or with B = 2 axial symmetry which were already discovered in the EYM theory. 

Finally, the numerical method employed to integrate the differential equations is based on the fourth-order Runge- 
Kutta method with a grid size 0.05. 


* Electronic address; norikoshiiki@mail.goo.ne.jp 
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FIG 



FIG. 1: The profile function / as a function of the radial coordinate x for |A| = 0.01, 0.1, 0.2 with a = 0.1 fixed. 



2: The solid line shows the a dependence of the maximum value of |A| above which there exists no soliton solution. 



FIG. 3: The dependence of the ADM mass M on the cosmological constant |A| for a = 0.06, 0.08 and 0.1. 










FIG. 4: The profile function / as a function of the radial coordinate x for |A| = 0.0, 0.5,1.0 with Xh = 0.1 and a = 0.02 fixed. 
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FIG. 5: The horizon radius Xh as a function of the black hole skyrmion mass Mhh for |A| = 0.0,1.0, 2.0 with a — 0.02 fixed. 
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FIG. 6: Shooting parameter fh as a function of |A| for a = 0.0, 0.02, 0.04. 
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FIG. 7: The solid line shows the a dependence of the maximum value of |A| above which there exists no black hole solution. 
The dotted line shows the a dependence of the value of |A| above which the lower branch solution change its stability to become 
stable. 



FIG. 8: Regular : The dependence of the ground state energy Wq on the cosmological constant [A| and the coupling constant 

a. 



FIG. 9: Black Hole : The dependence of the ground state energy Wq on the cosmological constant |A| and the coupling constant 
a with Xh = 0.1 fixed. 










